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ABSTRACT 

XMM-Newton observed the accreting millisecond pulsar SAX J1808. 4-3658 dur- 
ing its 2008 outburst. We present timing and spectral analyses of this observation, 
in particular the first pulse profile study below 2 keV, and the high-resolution spec- 
tral analysis of this source during the outburst. Combined spectral and pulse profile 
analyses suggest the presence of a strong unpulsed source below 2 keV that strongly 
reduces the pulsed fraction and a hard pulsed component that generates markedly 
double peaked profiles at higher energies. We also studied the high-resolution grating 
spectrum of SAX J1808. 4-3658, and found several absorption edges and Oxygen ab- 
sorption lines with whom we infer, in a model independent way, the interstellar column 
densities of several elements in the direction of SAX J1808. 4-3658. 
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1 INTRODUCTION 

The accreting millisecond X-ray pulsar (AMXP) SAX 
J1808.4-3658 (J1808 from now on) was the first X-ray bi- 
nary found to pulsate in the millisecond range, wi th a 
spin period of 2.5 ms jWiinands fc van der Klig|l998h ). It 
has been observed in outburst 6 times, roughly every 2.5 
years since 1996. During the outburst, the magnetic field is 
thought to channel part of the disc material onto the neu- 
tron star magnetic poles. The radiation emitted from the 
impact region (hot spot) and/or a slab of shocked mate- 
rial above it is then modulated at the neutron star spin 
period. This radiation is observed as pulsed emission that 
adds to the unpulsed emission coming from the accretion 
disc. A possible comptonizing medium surrounding the im- 
pact regi on can upscatter part of th e radiation to higher 
energies (|Poutanen fc Gierliriskill2003l ). The pulsations and 
the X-ray spectrum were observed during previous out- 
bursts by RXTE and a first st udy of the 1998 outburst was 
performed using t hose data ("Pouta nen fc GierlinskH |2003| . 
Ilbragimov fc Pouta ncn 2008). J 1808 was never observed be- 
low 2 keV during an outburst (except in th e 2000 and 2005 
out bursts at very low lum inosity levels, see IWiinandsllioOsI 
and lCampana et al.ll2008l ). This energy range is important 
to understand the pulse formation mechanism, because it is 
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here that both the hot spot and the accretion disc thermal 
emissions are expected to peak. Also, many absorption lines 
from the interstellar medium are expected in this energy 
range. These lines are important to definitively determine 
the interstellar column density toward the source. Broad- 
band spec tral analyses o f this XM M observation were re - 
ported bv iPapitto et all (|2009l ') and ICackett et al.l l|2009ll . 
who both focused on the study of the iron line emission at 
6-7 keV. Here we present the first simultaneous spectral and 
timing analysis of the pulsations of J1808 as observed with 
XMM-Newton during the 2008 outburst, with a particular 
attention to the lower energy range (< 2 keV). 



2 X-RAY OBSERVATION 

J1808 was observed in outburst with XMM-Newton on 2008 
October 1st (MJD 54740), for 63 ks of on-source exposure 
time. At the time of this observation, J1808 was at the 
beginning of the exponential decay stage of the outburst, 
with a relatively high flux level (see Hartman et al. 2009 
for a description of the over all outburst Ughtcu rve). The 
XMM-Newton Observatory (|jansen et all 120011 ') includes 
three 1500 cm'^ X-ray telescopes with the European Photon 
Imaging Camer a (EPIC), and a Reflec ting Grating Spec- 
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Energy [keV] Pulse phase [cycle] 

Figure 1. Left panel: fractional amplitude for the pulse profiles in 9 energy bands. The circles and triangles refer to the fundamental 
and second harmonic respectively. Vertical bars indicate 1 a errors in amplitude and horizontal bars the energy range. The black squares 
are the fractional amplitudes of the global profile, with the rms of the fundamental and second harmonic added in quadrature. The 
pulsed fraction is much lower below 2 keV, reaching a minimum value of ^ 0.4% rms in the 0.3-1 keV energy range. Right panel; 
pulse profiles in four different energy bands (blue asterisks 0.3-2 keV, red squares: 2—3.5 keV, green circles: 3.5-4.5 keV, black triangles 
4.5—12 keV). The profiles are double peaked, with the strength of the second harmonic increasing at higher energies. Two cycles are 
plotted, with profiles normalized to the maximum intensity and shifted by arbitrary amounts. The bottom profile corresponds to the low 
energies (0.3-2 keV) and its fractional amplitude is considerably reduced (although still significantly detected). 



work. The EPIC ca mera is compose d of two MPS CCD s 
l|Turner et al.l I2001I ) and a pn CCD (|Struder et al.ll2001^ . 
Each EPIC camera has a fixed, mode dependent frame read- 
out frequency, producing event lists in the 0.1-12 keV energy 
range. The RGS, is composed of a double array of gratings 
and produces high resolution spectra in the 0.33 to 2.5 keV 
energy range. 



Data have been processed using SAS version 8.0.1, and 
we have employed the most recent calibration files (CCF) 
available at the time the reduction was performed (February 
2009). Standard data screening criteria were applied in the 
extraction of scientific products. After removing solar fiares 
and telemetry dropouts the net pn exposure time is 41 ks. 
The central CCD of MOSl was operated in full frame 
mode with thin filters, and is heavily piled-up. For this 
reason we do not consider the MOSl data any further. The 
M0S2 was operated in timing window mode, and is also ex- 
cluded from our analysis since the 1.5 ms time resolution is 
insufficient to study the pulsations, and its spectral capabili- 
ties for very bright sources are not as much calibrated as the 
pn camera. The pn camera was operating in timing mode 
(with a thin filter), in order to reduce pile-up and allow 
the high precision timing analysis required for an accreting 
millisecond X-ray pulsar (AMXP) . We extracted the source 
photons from the pn with RAWX coordinates 26-49. The 
background is obtained from a region of the same size, at 
RAWX 2-25. Only photons with PATTERN^ 4 were used. 
The extracted spectrum was rebinned before fitting to ob- 
tain at least 100 counts per bin and the pn energy resolution 
was not oversampled by more than a factor three. We also 
extracted first and second order RGSl and RGS2 spectra, 
using the standard procedure reported in the XMM-Newton 
analysis manual, and again we ensured to have at least 100 
counts per spectral bin in any RGS spectrum. 



3 TIMING ANALYSIS 

We have first corrected the event times to the barycen- 
ter of the Solar System (us ing the SAS tool bar ycen, and 
the optical position given in lHartman et al]l2008l) and then 
we applied the 2008 o utburst timing solution published in 
iHartman et alj (|2009l ) in order to predict the phases of 
each photon detecte d and reconstruct the pulse profiles (see 
IPatruno et al] |2009| for a detailed explanation of the tim- 
ing technique). For the timing analysis we use events in the 
0.3-12 keV energy range. The pulse profiles are built by fold- 
ing data segments of length ~ 3500 s for the pulse phase 
analysis. This length is chosen to guarantee sufficiently high 
signal-to-noise profiles even if the pulsed fractions are small. 
The pulses are then decomposed by fitting two harmonics 
with frequency fixed at the pulse frequency (fundamental, 
v) and twice the pulse frequency (second harmonic, 21/) plus 
a constant representing the non-pulsed emission. 

We did not attempt to calculate a new timing solu- 
tion since the precision of the solution achievable with the 
short observation baseline of XMM is at least an order of 
magn itude lower than wha t was obtained with the RXTE 
data l|Hartman et al.ll2009l ). The short baseline of the ob- 
servation is also insufficient to model the timing noise that 
af fects the pulse phases and th at was extensively discussed 
in iHartman et al.l l|2008l . |2009| ) for J1808. If timing noise is 
present, a systematic error is introduced in the determina- 
tion of the pulse phases and spin frequency jPatruno et al.l 
12009.1 . 

Therefore we decided to subtract the solution reported 
in lHartman et all (|2009l ') and obtain the phase residuals with 
respect to that constant pulse frequency plus Keplerian cir- 
cular orbit model. The pulse phase residuals of the funda- 
mental drift by ~ 0.1 cycles during the observation. We also 
found a correlation between these pulse phase residuals and 
the 0.3-12 keV X-ray flux. We fitted the data with a linear 
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relation that gives a 11.4 for 10 degrees of freedom, 

and a slope of (9.9 ± 1.2) x 10~^ cycle/ct/s (we quote the la 
error) . 

The second harmonic is significantly detected in only 
ss 1/3 of the profiles, where a detection is defined as a ratio 
between the pulse amplitude and its statistical error larger 
than 3. When not detected, the second harmonic fractional 
amplitude upper limit was between 0.4 and 0.6% rms at the 
98% confidence level. No correlation between pulse phases 
and fiux was found for the second harmonic. 

To increase the signal to noise and calculate the har- 
monic content of the pulsations in the whole energy band, 
we folded the entire 41 ks data into one single pulse profile. 
The fractional rms amplitude of fundamental and second 
harmonic in the 0.3-12 keV energy band is 0.98(2)% rms 
and 0.33(3)% rms respectively {la uncertainties). 

We then repeated the procedure by dividing the obser- 
vations in 9 energy bands between 0.3 and 12 keV. The frac- 
tional amplitude of the pulse profile increases from 0.3 up 
to 3 keV, and then it remains constant within the errors up 
to 12 keV (Fig[T|. The fundamental tracks the behaviour of 
the total pulse profile. The second harmonic increases mono- 
tonically in the energy range considered. At energies above 
« 6keV the fractional amplitudes of the fundamental and 
second harmonic are comparable an d the overall pulse pro - 
file is double peaked (Fig. [l] see also lHartman et al.ll2009r ). 



4 SPECTRAL ANALYSIS 

We performed spectral analysis using the EPIC-pn spectrum 
(extracted as reported in § 2) in the 0.6-12 keV energy range, 
and the RGSl and RGS2 in the 0.4-1.8 keV energy range for 
the first order, and 0.7-1.8 keV for the second order. XSPEC 
version 11.3 was used for the spectral analysis. 

To account for relative flux calibration uncertainties be- 
tween different instruments, we fitted a multiplicative con- 
stant with the model, allowing for up to 10% relative calibra- 
tion flux scaling between EPIC-pn and RG^ The relative 
offsets between the instruments are found to be less than 
4% in all our flts. Furthermore, we included a 1.5% system- 
atic to the errors of each spectral bin (using the relative 
XSPEC tool) to take into account the calibration inaccuracies 
of each singl e instrument usecfl. We fir st used solar abun- 
danc es from Anders fc Grevessd (Il989l) and cross-sections 
from iBalucinska-Church fc McCammonI (|l992f ) for the pho- 
toelectric absorption. We tried an absorbed power-law plus 
a multi temperature disc and a single temperature black- 
body mod el (phabs* (dis k bb + bbody + powerlaw) ) as sug - 
gested in ICackett etlll (|2009| ) and IPapitto et ail l|2009r ). 
The — 5 unacceptable, mainly because of unmod- 
elled features in the data points between 0.3-2 keV and 6- 
7keV. 

The 6-7 keV energy range is where fluorescence lines 
of Fe are expected (iG corgg_fc Fabian 1991). Following 
ICackett et"all (|2009l ') and lPapitto et all ([20091) who claimed 
the detection of a broad iron Kq line in this energy range, we 
fltted this feature with a diskline model (see Tab. 1). We 




Figure 2. uF,, plot of the model used to fit the continuum. The 
power-law flux dominates over the other spectral components. 



^ http:/ /xmm2.esac.esa.int/docs/documents/CAL-TN-0052-5- 
O.ps.gz 

^ http:/ /xmm2.esac.esa.int/external/xmm_sw_cal/calib 



refer to ICackett et aP (|2009l ) and IPapitto et all l|2009l ) for 
discussion of this broad iron line. To model the 0.3-2 keV 
features, we flrst tried several photoelectric cross-sections, 
and different element abundances. The residuals are not very 
sensitive to the photoelectric cross-section parameters, while 
they strongly depend on the assumed abundances. The best 
flt model, however not yet satisfactorily, w as found using 
the lBalucinska-Church fc McCammonI l|l992h cross-sections, 
and the IWilms et al.l ()200GI 1 abundances. The stronger low 
energy features in the residuals were coming from absorption 
and emission features close to the oxygen K-edge (Ox-edge) 
at 0.543 keV. We decided to model only the continuum and 
the iron line as a first step, ignoring the data between 
0.5-0.6 keV. This energy range and the single features it con- 
tains were then investigated separately by using the RGS 
data (see § 14.11 and Table 1 for the continuum and Fe line 
spectral results). 

The two weak features at 1.8 keV and 2.2 keV are known 
to be due to the instrumental Si and An K edges, not yet 
perfectly calibrated, especially when dealing with timing 
mode observations (see the EPIC calibration report in the 
footnote). W e do not find evidenc e of a 0.871 keV O VII edge 
(reported in IPapitto et al.l |2009| ) , and no other significant 
edges are detected in the 0.8-1 keV range. 

After removing the 0-Si-Au edges (0.5-0.6 and 1.6- 
2.3 keV), we obtain xl = 1-41 (4051 dof). Given the high 
quality X-ray spectrum, the relatively high value of Xv = 
1.41 is very likely due to inter-calibration problems between 
the pn and the RGS spectra, and to the calibration uncer- 
tainties of each camera which emerge when observing bright 
sources. In fact, when using only the pn spectrum, we ob- 
tain a statistically acceptable fit with a = 0.95 for 235 
dof, with the same spectral parameters as reported in Tab. 1. 
The same model applied to the RGS alone gives a x^ = 1-1. 
Therefore we accept the xl = 1-41 and do not further com- 
plicate the spectral model. 



4.1 High-resolution spectroscopy 

We inferred the interstellar medium (ISM) abundances in 
a model independent way, by separately fitting the absorp- 
tion edges to the RGS data. We used only the RGSl and 
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RGS2 first order spectra, which are the best calibrated and 
have a higher number of counts. The continuum parame- 
ters were Icept fixed at the value reported in Tablet. We 
used vphabs that allows to set fixed abundance parameters 
with respect to the solar composition and to isolate the sin- 
gle absorption edges. The strongest features were observed 
around the Ox-edge {Ok; see Fig|3}. To model this edge we 
fixed at zero the oxygen abundance of the vphabs model, 
and fit only the data around the Ok with an edge model 
(note that the continuum is relatively constant in the small 
energy range around the edge itself, hence the modelling is 
independent on the broadband continuum model). The best 
fit gives Ok = 0.5421 ± 0.0003 keV and r = 0.66 ± 0.01, 
with Tx = Nx X a and g the photoelectric cross section. 
With the same method we fitted also the iron L [FeL), neon 
K (Nck), magnesium K [Mqk) and silicon K [SIk) edges 
that lie in the RGS band. The column density for each el- 
ement is reported in Tab[2] The absorption features close 
to Ok were fitted with Gaussian lines and were identified 
as Is — 2p atomic transitions of OI, Oil and OIII. The 
OIV line lies too close to an instrumental bad column to 
be detectable. The levels higher than IV are not detected, 
with a 98% upper limit of 0.08 eV on the equivalent width 
(EW). Similar features ha ve been observed in ot her X-ray 
bina ries like XB 1254-690 (iDiaz Trigo et al.ll2009l ') and Cyg 
X-2 (|Takei et al.ll2002l : ICostantini et al.ll2005l V 

From the single absorption edges we calculated the 
column density of each element (t x = Nx x a) assum- 
ing photoelectric cross sections from lGould fc June! l|l99l[) . 
The equivalent hydrogen column density is inferred from 
the best me as ured edge (Ok), by using abundances from 
IWilms et all l|2000l ). We derived Nh = (1.4 ± 0.2) x 
lO'^^ cm~^, consistent with the value derived in the direc- 
tion of J1808 from both HI and HII measurements (1.3 x 
10^^ cm~^ and 1.14 x 10 ^^ cm"^; iDickev fc LockmanI [l990l 
and lKalberla et ai]|2005l respectively). 

From the measurement of the EW of the oxygen I, 
II, and III Is — 2p transitions, using the method of the 
curve of growth, we have an independent measure of the 
relative oxygen abundances, using the formula EW\ = 
8.85 X 1Q~^'^ Nx^? fij , where A is the wavelength of t he line, 
and /ij the oscillator strength for the transition l|Spitzeij 
1 19781 ). We found an A'oxygon consistent with that inferred 
from the Oif-edge (see Tab. 2). 



5 DISCUSSION 

We observed for the first time the pulsations in the 0.3- 
2 keV energy range for J1808. The pulse fractional am- 
pHtudes sharply decrease below « 2keV during the 2008 
XMM-Newton observation. Our spectral analysis required 
a multi temperature blackbody at low energies that we iden- 
tify with the accretion disc emission, as well as the main 
source of unpulsed emission. We interpret the higher energy 
blackbody as the emission coming from the hot spot of the 
neutron star. The power law is also a pulsed component, 
since we measured ~ 2% rms pulse amplitudes up to 12 keV 
were the power law emission dominates. 

The empirical power-law model requires a disc con- 
tributing only « 30% of the power-law flux below 2 keV 
(Fig [2}. Therefore the disc is not the only responsible of the 
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Figure 3. High resolution RGSl spectrum of J1808 around the 
Ox-edge: OI, Oil and OIII absorption lines for the ls-2p transi- 
tion. The O IV line falls too close to the instrumental edge to be 
measured. 



Parameters 


DiskBB-(-BB+PL-|-Diskline 


Nh (10^2 cm-2) 


0.16 ±0.02 


Inner Disc kT (keV) 


0.20 ±0.01 


Disc flux (ergs~^cm~2) 


(1.89 ±0.21) X 10-10 


kTBB (keV) 


0.33 ±0.01 


BB radius (km) 


10.6 ±3.2 


BB flux (ergs~^cm~2) 


(1.22 ±0.15) X 10-1" 


Photon Index F 


2.11 ±0.01 


PL flux (erg s~^cm~'^) 


(1.53 ± 0.09) X 10-8 


Ep.e (keV) 


6.45 ±0.08 


EW (eV) 


97.7 ±31.4 


R/jv (km) 


20 ± 2 


RoiJT (km) 


193 ± 15 


Incl. (deg) 


> 44" 




-2.1 ±0.1 


Fe-line flux (ergs ^cm 2) 


(6.8 ±0.2) X 10-" 


Flux (ergs^^cm"'^) 


(1.85 ± 0.08) X 10-9 


Absorbed Flux (ergs~^cm~2) 


(1.51 ± 0.08) X 10-9 


xl (dof) 


1.41 (4051) 



Table 1. Spectral parameters for J1808 in outburst, from com- 
bined pn and all RGS data. Errors ar e at Icr confldence level. 
N^f is calculated with abundances from I Wilms et aL ] 1I2OOOI) . Un- 
absorbed fluxes are given in the 0.5-10 keV energy range. The 
blackbody radius is calculated assuming a distance of 3.5 kpc and 
a neutron star mass of 1.4Mq . a. The lower limit is quoted at 95% 
confldence level, b. /3 is the power law index of the emissivity 



steep decrease of pulse fractional amplitudes below 2 keV. 
However, an empirical power-law model is unphysical and a 
self consistent physical scenario would require a sharp de- 
crease of the power-law component below ft; 2 — 3 keV. In 
that case the unpulsed disc emission might dominate below 
2 keV and explain the sudden drop of the pulse amplitudes 
at those energies. A Comptonization model with a shocked 
plasma in a slab geometry is expect ed to cut off below fti 3 
keV a nd was already proposed by iGierliilski fc PoutanerJ 
(|2005l ) for the AMXP XTE J1751-305. The hot single tem- 
perature blackbody would then be produced by the pulsat- 
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Edge 


Energy fkeV) 


T 


a 






n ^A9i -1- n nnn'i 

U.O^^l ^ u.uuuo 


n fifi -t- n ni 




11 7-1-0 9 


FeL 


0.712 ±0.005 


0.08 ±0.02 


4.936 


1.6 ±0.4 


NcK 


0.865 ± 0.004 


0.10 ±0.02 


3.523 


2.8 ±0.5 


Mgx 


1.281 ±0.007 


0.06 ±0.01 


2.191 


2.7 ±0.4 


Six 


1.79 ±0.007 


0.14 ±0.01 


1.476 


9.4 ±0.6 


ls-2p 


Energy (keV) 


EW (eV) 


^oxygen 




OI 


0.527 ±0.001 


1.2 ±0.2 


9 ± 2 




Oil 


0.531 ±0.001 


1.2 ±0.2 


8 ± 2 




OIII 


0.537 ±0.001 


0.6 ±0.1 


11 ± 2 





Table 2. Interstellar medium edges and lines in the direction of 
SAX J1808. 4-3658. Photoelectric cross sections (tr) are in units of 
10~^^ cm^ (from Gould & Jung 1991), while the column densities 
Nx are in units of lO^^cm"^. 



ing emission of the hot spot, the multi temperature black- 
body by the unpulsed radiation of the accretion disc and the 
hard component by the pulsed comptonized radiation of the 
shock around the hot spot. 

The origin of the second harmonic can be related to a 
different pulse profile of the hard emission compared to the 
lower- energ y blackbody as was suggested by lGierliriski et al.l 
l|2002h and lPoutanen fc Gierlihskr(l2003f) for RXTE ob ser- 
vations of J1808. and by iGierlihski fc PoutanenI (j2005l ) for 
XMM and RXTE obseryations of XTE J1751-305. This 
might result from a different angular distribution of the 
Comptonized radiation, which is expected as it is produced 
in the optically thin accretion shock, but not at the stellar 
surface as the blackbody emission. 

The reason why the phase of the fundamental is cor- 
related with the X-ray flux while the second harmonic is 
not, can then be related with different formation processes 
for the fundamental and the second harmonic. If the hot 
spot contributes only to the fundamental frequency while 
the comptonization region contributes to both the funda- 
mental and the second harmonic, the fundamental pulse 
phase may track the hot spot position. The phase of the 
second harmonic instead will be affected by the comp- 
tonization process and might come from an extended re- 
gion around the hot spot, coyering a large area of the neu- 
tron star surface. The hot spot can instead come from a 
well defined region on the neutron star surface, and can 
moye according to the X-ray flux fluctuations thus produc- 
ing the pulse pha se wandering correlated with X-ray flux 
ijRomanoya et al.l [2004. Lamb ct al. 2008). 

We also measured the first model independent column 
densities of several elements in the line of sight of J 1808 
(Tab. 2). The most precise measurement comes from the oxy- 
gen column density, from which we could deriye the equiva- 
lent hydrogen colum n density Nh = 1.4 x 10^^ cm~^, assum- 
ing abundances from I Wilms et aDl2000l . This determination 
is particularly important: i) in the study of the X-ray emis- 
sion of this object during quiescence and cooling, where the 
uncertainty in the assumed Nh value could alter the reli- 
ability of the sour ce intrinsic luminosity (see for example 
iHeinke et aLll2009l for coo ling studies of SAX J1808.4-3658 
and lYakovlev et al ] |2005l for a discussion of the problem), 
and ii) for optical studies who can now rely on a more pre- 
cise determination of the extintion value toward this system. 



It is ineresting to compare the abundances we derived 
from the ed ge fitting with the expected value for the ISM (as 
reported bv lWilms et al. P2000Y We find that in the direction 
of J1808, the Ne/O abundance is ~0.23, slightly larger than 
in the ISM (~ 0.18), which maybe be poiting to a Ne-rich 
environment as observed in other low mass X-ray binaries 
(Juett et al. 2001). 
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